Identifying the viscoelastic properties of soft matter from the indentation response of a hard film-soft substrate system

The indentation technique is widely used in measuring the mechanical properties of soft matter at the microscale or nanoscale,but still faces challenges by these unique properties as well as the consequent strong surface adhesion, including the strong nonlinear effect, unclear judgment of the contact point, difficulties in estimating the contact area, and the risk of the indenter piercing the sample. Here we propose a two-step method to solve these problems: lay a hard film on a soft matter, and obtain the viscoelastic properties of this soft matter through the indentation response of this composite structure. We first establish a theoretical indentation model of the hard film-soft substrate system based on the theory of plates, elastic-viscoelastic correspondence principle and Boltzmann superposition principle. To verify the correctness of this method, we measure the mechanical properties of the methyl vinyl silicone rubber(MVSR) covered by a Cu nanofilm. Finally, we test the effectiveness and error sensitivity of this method with the finite element method(FEM). The results show that our method can accurately measure the mechanical properties of soft matter, while effectively circumventing the problems of the traditional indentation technique.     

A general modal frequency-domain vortex lattice method for aeroelastic analyses

The generalized aerodynamic force (GAF) matrix is derived for the Unsteady Vortex Lattice Method (UVLM) without the assumption of out-of-plane dynamics. As a result, the approach naturally includes in-plane motion and forces unlike the doublet lattice method (DLM). The derived UVLM GAF is therefore applicable to industry-standard techniques for aeroelastic stability analyses, such as the p–k method. In this work, the fluid–structure interpolation is performed with radial basis functions for surface interpolation. The generalized aerodynamic forces computed with the UVLM are verified against the DLM from NASTRAN on a simple flat plate configuration. The ability of the UVLM to include steady loads is verified with a T-tail flutter case and the results confirm the importance of including steady loads for T-tail flutter analysis. The modal frequency domain VLM therefore provides the same level of efficiency and accuracy than the DLM, but without the restrictions and with the ability to handle complex geometries. It is therefore a viable replacement to the DLM.     

基于局部加密纯无网格法非线性Cahn-Hilliard方程的模拟

为数值求解描述不同物质间相位分离现象的高阶非线性Cahn-Hilliard(C-H)方程,发展了一种基于局部加密纯无网格有限点集法(local refinement finite pointset method,LR-FPM).其构造过程为:1)将C-H方程中四阶导数降阶为两个二阶导数,连续应用基于Taylor展开和加权最小二乘法的FPM离散空间导数;2)对区域进行局部加密和采用五次样条核函数以提高数值精度;3)局部线性方程组求解中准确施加含高阶导数Neumann边值条件.随后,运用LR-FPM求解有解析解的一维/二维C-H方程,分析粒子均匀分布/非均匀分布以及局部粒子加密情况的误差和收敛阶,展示了LR-FPM较网格类算法在非均匀布点情况下的优点.最后,采用LR-FPM对无解析解的一维/二维C-H方程进行了数值预测,并与有限差分结果相比较.数值结果表明,LR-FPM方法具有较高的数值精度和收敛阶,比有限差分法更易数值实现,能够准确展现不同类型材料间相位分离非线性扩散现象随时间的演化过程.     

A Runge-Kutta–based WENO reconstruction on moving mesh for compressible fluids

In this paper, we construct a high-order moving mesh method based on total variation diminishing Runge-Kutta and weighted essential nonoscillatory reconstruction for compressible fluid system. Beginning with the integral form of fluid system, we get the semidiscrete system with an arbitrary mesh velocity. We use weighted essential nonoscillatory reconstruction to get the space accuracy on moving meshes, and the time accuracy is obtained by modified Runge-Kutta method; the mesh velocity is determined by moving mesh method. One- and two-dimensional numerical examples are presented to demonstrate the efficient and accurate performance of the scheme.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

A framework to perform high-order deconvolution for finite-volume method on simplicial meshes

In this article, a new framework to design high-order approximations in the context of node-centered finite volumes on simplicial meshes is proposed. The major novelty of this method is that it relies on very simple and compact differential operators, which is a critical point to achieve good performances in the High-performance computing context. This method is based on deconvolution between nodal and volume-average values, which can be conducted to any order. The interest of the new method is illustrated through three different applications: mesh-to-mesh interpolation, levelset curvature computation, and numerical scheme for convection. Higher order can also be achieved within the present framework by introducing high-rank tensors. Although these tensors feature much symmetries, their manipulation can quickly become an overwhelming task. For this reason and without loss of generality, the present articles are limited to third-order expansion. This method, although tightly connected to the k-exact schemes theory, does not rely on successive corrections: the high-order property is obtained in a single operation, which makes them more attractive in terms of performances.     

L1-type smoothness indicators based weighted essentially nonoscillatory scheme for Hamilton-Jacobi equations

This article presents an improved fifth-order finite difference weighted essentially nonoscillatory (WENO) scheme to solve Hamilton-Jacobi equations. A new type of nonlinear weights is introduced with the construction of local smoothness indicators on each local stencil that are measured with the help of generalized undivided differences in L¹-norm. A novel global smoothness measurement is also constructed with the help of local measurements from its linear combination. Numerical experiments are conducted in one- and two-dimensions to demonstrate the performance enhancement, resolution power, numerical accuracy for the proposed scheme, and compared it with the classical WENO scheme.     

核电结构土-结相互作用分析分区混合计算方法

土-结构相互作用分析是核电结构进行抗震设计和安全评估的重要环节.在核电结构的土-结相互作用分析中,阻尼和非线性是影响结构反应的重要因素. 若采用频域分析,可以方便考虑阻尼,但需通过等效线性化来考虑非线性,不适合于强震作用下的土体非线性.若采用时域分析的逐步积分方法,适合于考虑非线性,但材料阻尼一般采用瑞利阻尼模型,除了紧靠指定阻尼比的少数几个振型外,其他振型的反应将受到瑞利阻尼模型所确定的大阻尼所抑制,造成地震反应与真实情形有较大差异.若采用时域分析的模态叠加法,可合理计入阻尼效应,但模态叠加法不能考虑非线性.因此,如何合理考虑阻尼和非线性是核电结构土-结相互作用分析需要关注的问题.基于此,本文提出一种模态叠加和时步积分结合的土-结相互作用分区算法.其中,出于安全性考虑,地震作用下核电主体结构一般不允许进入非线性,因此结构可采用模态叠加方法,以便合理考虑结构阻尼;土体和基础采用显式时步积分法,可考虑土体非线性;通过人工边界条件考虑无限域的影响 (辐射阻尼).通过简单算例对该方法进行了验证,并用于CAP1400核电结构的土-结相互作用分析中,对比分析了采用模态阻尼和瑞利阻尼时核电结构和场地反应的差异,结果表明结构阻尼模型对场地的反应影响不大,但对结构反应影响明显,在实际工程中应合理选取阻尼模型.      

基于内聚力模型的快速荷载下CFRP-混凝土界面剪切剥离模拟

鉴于纤维增强复合材料(Fiber Reinforced Polymer, FRP)已在混凝土结构加固中普遍应用, 能有效提高结构的承载和变形能力. 为揭示动载作用下FRP-混凝土界面的剥离机理, 本文基于ABAQUS平台, 采用内聚力模型模拟CFRP-混凝土界面层, 实现了快速荷载下CFRP-混凝土界面剥离的高效模拟. 结果表明: CFRP表面应变在加载过程中由加载端向自由端方向传递, 随着加载速率的提高, 界面承载能力也随之提高; 界面峰值剪应力也存在显著的应变率效应; 模拟结果与试验结果基本吻合, 说明了新方法的有效性.